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ABSTRACT 


System  identification  techniques  were  used  to  estimate  muscle  forces 
during  horizontal  saccadic  eye  movements  in  order  to  better  understand  the 
neuronal  control  strategy.  The  lateral  and  medial  rectus  muscle  of  each  eye 
was  modeled  as  a  parallel  combination  of  an  active  state  tension  generator 
with  a  viscosity  and  elastic  element,  connected  to  a  series  elastic  element. 
The  eyeball  was  modeled  as  a  sphere  connected  to  a  viscosity  and  elastic 
element.  The  predictions  of  the  model  were  shown  to  be  in  good  agreement  with 
the  data.  The  results  of  extensive  analysis  did  not  support  the  existence  of 
a  postulated  continuous-time  external  feedback  control  mechanism.  Analysis  of 
the  data,  however,  did  support  a  time  optimal  control  strategy,  a  strategy 
which  directs  the  eyeball  to  its  destination  in  minimum  time  for  saccades  of 
all  sizes.  1 


I 


Saooadio  eye  movements  are  among  the  fastest  voluntary  muscle  movements 


the  human  body  Is  capable  of  producing  and  are  characterized  by  a  rapid  shift 
of  gaze  from  one  point  of  fixation  to  another.  Although  the  purpose  for  such 
an  eye  movement  is  obvious/ !  that  is,  to  quickly  redirect  the  eyeball  to  the 

target^  the  neuronal  control  strategy  is  not. -  Por  instance,  does  the  word 

/ 

"quickly"  in  the  previous  sentence  imply  the  most  rapid  movement  possible,  or 
simply  a  fast  as  opposed  to  a  slow  movement?  7  The  present  investigation 
utilizes  system  identification  techniques  to  estimate  muscle  forces  during 
horizontal  saccadic  eye  movements  in  order  to  better  understand  the  neuronal 
control  strategy.  . 

Until  quite  recently,  models  of  the  saccadic  eye  movement  system  involved 
a  ballistic  or  preprogrammed  control  to  desired  eye  position  based  on  retinal 
error  alone.  (1)  -  (4)  Today,  an  increasing  number  of  authors  are  putting 
forth  the  idea  that  goal-directed  saccades  are  controlled  by  a  local  feedback 
loop  that  continulously  drives  the  eye  to  the  desired  position.  This  hypo¬ 
thesis,  first  presented  by  Vossius  in  I960,  did  not  start  to  gain  acceptance 
until  1975  when  Robinson  re-examined  it.  (5) (6)  Robinson  suggested  that 
saccades  originated  from  neural  commands,  which  encode  saccade  velocity  and 
duration,  to  the  pulse  generator  that  specifies  the  desired  position  of  the 
eye  rather  than  the  preprogrammed  distance  the  eye  must  be  moved.  The  image 
of  the  internal  representation  of  the  current  eye  position  is  subtracted  from 
the  internal  representation  of  the  desired  position  to  create  an  error  signal 
that  completes  the  local  feedback  loop  to  generate  the  neural  control  signal. 


This  neural  pulse  continiously  drives  the  eye  until  the  internal 
representation  of  the  error  signal  is  zero. 

One  of  the  basic  tasks  of  this  project  is  to  investigate  whether  an 
external*  continuous-time  sensor  feedback  control  nechanism  operates  during  a 
saocadic  eye  movement.  Studies  by  other  investigators  indicate  that  the 
profile  of  a  saccade  can  be  altered  up  to  50  msec  after  target  movement*  and 
that  visual  information  is  processed  continuously  during  this  internal.  (7) 
(8)  Their  results  indicate  that  50  msec  after  target  movement*  the  saccade 
profile  can  not  be  altered  by  altering  the  target  location.  The  present 
experiment  attempts  to  verify  their  results  by  artifically  maintaining  a 
constant  error  throughout  the  saccade  by  moving  the  target  the  same  distance 
the  eyeball  moves.  This  method  of  artifically  openning  the  feedback  loop  was 
used  by  Collewijn  and  Van  der  Mark  in  their  study  of  the  slow  phase  of 
optokinetic  nystagmus.  (9) 

A  second  task  in  this  project  involves  the  investigation  of  a  horizonatal 
saccadic  eye  model  presented  by  Enderle.  (10)  Estimated  values  of  the  poles 
of  the  transfer  function  showed  an  increasing  trend  as  target  displacement 
increased.  As  stated,  a  number  of  reasons  may  account  for  this,  a 
significantly  nonlinear  model  and/or  an  inaccurate  model,  or  poor  initial 
estimates  for  the  parameter  values.  It  is  essential  that  the  model  for 
horizontal  saccadic  eye  movements  be  accurate  since  it  is  to  be  used  in 
sensitively  assessing  whether  an  external  feedback  control  mechanism  operates 


during  saccade s. 


II.  OBJECTIVES 


The  main  objective  of  this  project  ie  to  better  understand  the  mechanisms 
controlling  saccadic  eye  movements.  The  specific  objectives  were: 

1.  Update  and  improve  the  model  for  horizontal  saccadic  eye 
aoveaents. 

2.  Estiaate  model  paraaeters  using  system  identification 
techniques. 

3.  Investigate  the  hypothetical  operation  of  an  external 
feedback  control  aechanisa. 

III.  MODEL  OP  OCULOMOTOR  PLANT 

In  the  development  of  the  model  for  horizontal  saccadic  eye  aoveaents  by 
Enderle.  it  vas  assumed  in  equation  2  that  the  time  derivative  of  the 
antagonist  tension  vas  negligibly  small.  This  assumption  vas  a  major  factor 
in  the  resulting  functional  relationship  betveen  target  displacement  and  poles 
of  the  transfer  function.  Another  faotor  effecting  this  relationship  vas  the 
initial  parameter  estimates,  vhich  vill  be  discussed  in  the  next  section. 
Using  Laplace  variable  analysis  about  the  operating  point  or  initial  eye 
position,  it  is  possible  to  derive  the  folloving  nev  and  more  accurate  model 
of  horizontal  saccadic  eye  movements  appropriate  for  the  system 
identification  technique  (See  Appendix  A1  for  details  of  the  derivation  and 
paraaeter  description).  (11) 

‘'WW'WmW  -  •  *  V  ♦  p25  *  P1« +  V-  (1> 

This  model  is  a  modification  and  correction  of  the  linear  homeomorphic  aodel 
by  Bahill  et.  al.  (4)  The  implementation  of  this  model  in  the  parameter 
estimation  routine  resulted  in  an  excellent  match  betveen  model  prediction  and 
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the  data.  No  apparent  trends  between  the  poles  of  the  transfer  function  and 
target  displacement  were  observed. 

Presented  in  Fig.  1  is  a  block  diagram  of  the  saccadic  eye  movement 
system  which  functions  as  one  element  within  the  overall  system  controlling 
all  oculomotor  movements.  The  saccadic  eye  movement  system  consists  of  three 
components;  the  controller,  that  is.  the  agonist  and  antagonist  active  state 
tension  generators,  the  oculomotor  plant,  and  the  feedback  element  H.  The 
oculomotor  plant,  indicated  by  the  section  within  the  dashed  lines  in  Fig.  2, 
is  based  on  the  Laplace  transform  of  equation  1.  Note  that  under  an  assumed 
ballistic  control  of  saccadic  eye  movements,  time  delays  in  the  feedback 
element  cause  this  system  to  operate  in  an  open  loop  mode.  If,  in  fact,  an 
external  continuous- time  feedback . control  mechanism  exists,  then  openning  the 
feedback  loop  causes  the  system  to  overreact  because  of  the  constant  error 
present  at  the  input. 

Since  no  assumptions  concerning  the  saccadic  active  state  tensions  are 
required  in  the  derivation  of  the  modified  linear  homeomorphic  model,  any  type 
of  control  signal  may  be  utilized.  For  example,  investigators  may  utilize 
this  model  when  examining  alternatives  to  the  pulse-step  muscle  force.  If  a 
continuous-time  external  feedback  control  mechanism  exists,  then  a  new 
waveform  could  be  postulated  for  the  active  state  tension.  This  waveform 
could  then  be  revised  based  on  the  predicted  response  of  the  modified  linear 


homeomorphic  model 


The  mechanical  components  of  the  model  are  conncected  such  that  the  rate 
of  change  of  the  mnscle  forces  together  with  the  actual  muscle  forces  drive 
the  eyeball  to  its  final  position.  Consider  the  lov-p&ss  filtered  pulse-step 
active  state  tension  signals  presented  in  Pig.  2.  During  the  intial  pulse 
phase  of  the  trajectory,  the  agonist  active  state  tension  drives  the  eyeball 
and  the  antagonist  active  state  tension  restrains  the  eyeball.  Since  the 
antagonist  active  state  tension  falls  to  zero  during  this  interval,  then  the 


rate  of  change  antagonist  active  state  tension,  F^j,,  is  negative  and  as 
indicated  by  the  equation  1,  acts  to  drive  the  eyeball  along  with  the  agonist 
force,  the  agonist  rate  of  change  force,  F^,.  As  a  result,  a  greater 
force  propels  the  eyeball  than  if  only  the  agonist  and  antagonist  active  state 
tensions  are  used  as  shown  in  the  plant  input  curves  of  Fig.  2.  This  result 
is  in  agreement  with  the  tension  data  recorded  for  horizontal  saccadic  eye 
movements  (12) (13).  After  the  pulse  phase  of  the  trajectory,  the  agonist 
active  state  tension  falls  to  a  steady  state  value  which  results  in  a  negative 
rate  of  change  agonist  active  state  tension,  F^,.  Thus,  according  to  equation 


1,  the  agonist  rate  of  change  force,  F^,,  acts  to  restrain  the  eyeball  with 


the  antagonist  force,  F^j.,  and  the  antagonist  rate  of  change  force,  F^j.. 


As 


shown  for  the  plant  input  curves  in  Fig.  1,  the  resultant  forces  act  like  a 
dynamic  brake  which  is  reduced  to  zero  as  the  eyeball  reaches  its  destination. 


IV.  SYSTEM  IDENTIFICATION  TECHNIQUE 


Paramter  estimates  for  the  model  of  the  oculomotor  system  presented  in 
the  previous  section  are  found  using  the  system  identification  technique  as 
described  by  Enderle.  (10) (14)  The  system  identification  technique  is  a 
frequency  response  method.  Ideally,  the  transfer  function  is  equal  to  the 


Fig.  2  Transformation  of  the  active  state  tensions  by  the  numerator 
terms  of  the  ocnlomotor  plant.  One  diagram  shows  the  assumed 
low-pass  filtered  pulse-step  waveform  of  the  agonist  and 
antagonist  active  state  tensions  initiated  at  100  msec.  The 
antagonist  pulse  circuiscribes  the  agonist  pulse  by  3  msec  on 
each  end.  Parameter  values  are  from  Bahill  et.  al.  for  simul- 
ting  a  10  degree  horizontal  saccadic  eye  movement  (4).  The 
other  diagram  shows  a  parital  and  full  transformation  of 
the  active  state  tensions  by  the  numerator  terms  of  the  ocul¬ 
omotor  plant. 


ratio  of  the  Fourier  transform  of  the  ontpnt  to  the  system  input.  For  the 
oculomotor  system,  the  transfer  function  is  calculated  from  the  fast  eye 
response  to  a  step  target  displacement.  Because  of  unequal  time  delays  and 
eye  displacements  from  saccade  to  saecade  for  the  same  target  displacement,  it 
is  not  possible  to  use  averaging  techniques  to  reduce  the  effects  of 
measurement  noise.  Fortunately,  the  measurement  noise  is  small  relative  to 
the  input  and  output  signals,  so  the  second  method  described  by  Otnes  and 
Enochson  is  used.  (15)  The  fast  eye  response  measurements  are  first  filtered 
using  a  Butterworth  low-pass  filter  with  a  half-power  point  at  125  hertz. 
Transforming  the  filtered  measurements  directly  by  the  Fast  Fourier  algorithm 
resulted  in  distortion  due  to  truncation,  since  the  signal  did  not  go  to  zero 
at  steady  state.  This  is  circumvented  by  subtracting  the  steady  state  value 
from  each  sample,  passing  this  signal  through  a  Kaiser  window,  packing  with 
zero's,  and  Fast  Fourier  transforming  the  modified  signal.  The  Fourier 
transform  of  the  fast  eye  response  is  equal  to  the  Fast  Fourier  transform  of 
the  modified  signal  plus  the  Fourier  transform  of  the  nnit  step  with  amplitude 
equal  to  the  steady  state  value.  The  input  signal  is  the  Fourier  transform  of 
the  unit  step  function  with  amplitude  equal  to  the  steady  state  value. 
Parameter  estimates  for  the  oculomotor  model  are  calculated  using  the 
conjugate  gradient  search  program  similar  to  Seidel's,  which  minimizes  the 
integral  of  the  absolute  value  of  the  error  squared  between  the  model  and  the 
data.  (Id) 

Data  were  collected  from  subjects  seated  before  a  target  display  of  seven 
small  red  light  emitting  diodes  (LED),  each  separated  by  five  degrees.  The 
subject's  head  was  restrained  by  a  bit-bar.  The  subject  was  instructed  to 


follow  the  "Jumping"  target  which  moved  from  the  center  position  to  any  of  the 
other  LED's  and  then  returned  to  the  center  position.  The  order  of  the  target 
positions  as  well  as  the  time  interval  between  displacements  were  randomized. 
Data  were  only  recorded  for  the  initial  displacement  from  the  center  position. 

Horizontal  eye  movements  were  recorded  from  each  eye  using  an  infrared 
signal  reflected  from  the  anterior  surface  of  the  cornea-scleral  interface 
with  instrumentation  described  by  Engelken  et.  al.  (17).  Signals  for  both 
eyes  tracking  were  digitized  using  the  analog/digital  converter  of  the  DECLAB 
PDP  11/34  computer  and  stored  in  disk  memory.  These  signals  were  sampled  at  a 
rate  of  1000  samples  per  second  for  one-half  second. 

Presented  in  Fig.  1  is  a  block  diagram  of  the  oculomotor  system,  which  is 
suitable  for  the  case  of  a  ballistic  controller,  or  as  Indicated,  with  an 
external  continuous-time  feedback  control  configuration.  Under  the  assumption 
of  ballistic  control,  the  system  operates  in  an  open  loop  mode  due  to  time 
delays  in  the  feedback  element.  Under  the  assumption  of  an  external 
continlous-time  controller  during  experiments  in  which  the  feedback  loop  is 
artifically  openned,  the  system  also  operates  in  an  open  loop  mode.  Thus, 
the  analysis  of  the  data  is  performed  with  an  open  loop  system  regardless  of 
the  type  of  control  assumed  operating.  However,  the  signals  describing  the 
agonist  and  antagonist  force  generators  are  quite  different  dependent  on  the 
type  of  controller. 

Under  an  assumed  ballistic  control  mechanism,  the  agonist  and  antagonist 
active  state  tensions  are  described  with  the  following  low-pass  filtered 


pulse-step  waveforms 
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initial  magnitude  of  the  agoniat  active  atate  tension, 

pnlse  magnitude  of  the  agonist  active  state  tension, 

step  magnitude  of  the  agonist  active  state  tension, 

initial  magnitude  of  the  antagonist  active  state  tension, 

step  magitnde  of  the  antagonist  active  state  tension, 

latent  period,  the  time  interval  between  the  target  movement  and 

the  start  of  the  eye  movement, 

duration  of  the  agonist  pnlse  active  state  tension, 
activation  time  constant. 


dr 


=  deactivation  time  constant. 

The  steady  state  antagonist  active  state  tension  terms,  F 
removed  from  the  analysis  since 
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according  to  equations  7  and  8  in  Appendix  Al. 

Physiological  data  was  used  to  estimate  all  of  the  initial  parameter 


estimates.  (4) (13) (18)  Specifically 


K  «  66.4  Nt/m 


(Robinson) 


Kjj  -  32.  Nt/m 

(Robinson) 

Rse  -  125  Nt/m 

(Collins) 

“  3.5  Nt  see/m 

(Bahill) 

■  1.66  Nt  sec/m 

(Bahill) 

B  ■  3.1  Nt  sec/m 

(Collins) 

J  -  2.2xl0~3  Nt  sec2/m 

(Robinson) 

The  value  of  K  is  determined  from  the  steady  state  agonist  and  antagonist 
aroscle  tensions 


where 


and 


AG 


-  T 


ANT 


Kx 


»  17  +  0  gm  tension 


(2) 


T^j  *  17  -  0.3  0  gm  tension  (3) 
from  Fig.  1  of  Robinson  et.  al.  (18)  The  valne  of  the  muscle  viscosity  terms 
reported  by  Bahill  et.  al.  are  modified  according  to  their  incorrectly  derived 
result.  (4) 


BAG  "  KST  ®AG  /  KSE 
BANT  “  kst  bant  /kse 


Substituting  these  values  into  the  expressions  for  6,  PQ,  P^,  Pj  and  Pj  found 
in  Appendix  Al,  yields 


6  -  5.0938109  x  107 


P3  -  1.548  x  1<T 

P2  =  3.4455  x  105 

Pj  -  1.9724  x  107 

PQ  -  2.2631  x  10® 

The  initial  extimates  for  the  four  real  eigenvalues  as  determined  from  the 
eharaeteristic  equations  for  this  system  are  -15.27.  -66.1.  -173.36,  and 
-1293.28. 


Now,  the  steady  state  active  state  tensions,  P^GS'  FANTS*  FANTO*  *re 

computed  froa 


fag  -  F  tag  +  hru  -  *P1» 

SB 

and 

FANT  “  Is1  TANT  "  *LT(x  +  XP1)# 

SB 

derived  from  Fig.  1  of  Appendix  Al.  Substituting  equations  2  and  3  into  the 
above  equations  with  Xp^  equal  to  3.14  ma  (assumed  from  Robinson)  yields 

FA(J  -  .1089  +  .01850  NT 

and 

FANT  "  .1089-9.838  x  10_30  NT. 

Note  that  when  FAfJT  at  steady  state  is  greeted  than  11°,  then  F ANT  equals  0. 
The  initial  estimate  of  the  pulse  magnitude  of  the  active  state  tension,  fagp* 
is  0.9806  (a  maximum  innernation  signal  of  100  gms) .  (13) 

The  initial  estimate  of  the  duration  of  the  agonist  pulse  active  state 
tension  is  estimated  directly  from  the  data.  The  maximum  velocity  is  computed 
using  a  two— point  central  difference  derivate  method.  (19)  The  initial 


estimate  of  t  is  than  the  tiae  interval  between  the  start  of  the  saocade  and 
P 

the  tiae  at  maximua  velocity. 

The  initial  estinate  of  the  deactivation  tiae  constant  is  also 
estinsted  directly  froa  the  data  aa  the  tiae  interval  between  aazisraa  velooity 
and  the  end  of  the  saecade  divided  by  4.  The  activation  tiae  constant  equals 

2  asec.  (20) 

Great  care  was  exercised  in  evaluating  initial  paraaeter  estimates 
since  large  differences  froa  the  true  values  conld  cause  the  estiaation 
rontine  to  converge  to  snboptiaal  and  nonphysiologically  consistent  results. 
The  precision  of  the  paraaeter  estimation  rontine  results  presented  in  Pigures 

3  and  4  for  a  10°  target  movement  are  typical  for  the  three  subjects  tested. 
There  was  no  significant  trend  in  the  values  of  the  eigenvalues  with  target 
displacement  and  the  estimates  correlated  well  with  their  physiological  data 
derived  values.  Listed  in  Table  1  are  the  final  estimated  values  for  the 
eigenvalues  for  target  displacements  of  5.  10,  and  15  degrees. 

V.  EXTERNAL  FEEDBACK  CONTROL  MECHANISM 

Experiments  were  conducted  in  which  the  hypothesized  continuous-time 
external  feedback  control  mechanism  was  artifically  openned.  Data  was 
collected  on  three  subjects  in  two  stages.  First  stage  data  consisted  of 
natural  saecades,  elicited  from  target  displacements  of  -15  to  15  degrees. 
The  second  stage  consisted  of  data  collected  by  artifically  maintaining  a 
constant  error  of  5°  throughout  the  saccade  by  moving  the  target  (under 
eoaputer  control)  the  same  distance  the  eyeball  moved.  This  type  of  movement 
will  be  referred  to  as  the  modified  target  trajectory.  The  experiment  started 
by  the  target  moving  ±  5  degrees.  The  subject  responded  by  following  the 


TABLE  1 


Final  parameter  eatiaates  for  the  four  real  eigenvalues  for 
target  displacements  of  5,  10,  and  IS  degrees  left  of  center  position 

TAKGST  DISPLACEMENTS 


5° 

o 

o 

ft 

15° 

T  EXE 

Eigenvalue  1 

-14.4 

-13.9 

-13.7 

Eigenvalue  2 

-83.8 

-84.4 

-93.6 

Eigenvalue  3 

-190.9 

-174.6 

-157.4 

Eigenvalue  4 

-1298.7 

-1283.2 

-1259.4 

!HT  BXB 

Eigenvalue  1 

-14.8 

-14.4 

-12.35 

Eigenvalue  2 

-80.8 

-84.9 

-88.2 

Eigenvalue  3 

-185.8 

-165.4 

-160.1 

Eigenvalue  4 


-1302.2 


-1261.6 


-1223.7 


j naping  target.  Instantaneously,  et  the  eyeball  moved,  the  target  moved  that 
same  distance.  Significant  differences  in  saccade  profile  between  the  two 
stages  wosld  indicate  the  possible  existence  of  an  external  feedback 
mechanism.  Data  was  first  analyzed  using  the  two-point  central  difference 
method,  and  then  nsing  the  system  identification  technique. 

The  results  of  the  analysis  of  the  data  on  the  three  subjects  nsing  the 
two-point  central  difference  method  is  contained  in  Appendix  A2.  (19)  The 
results  of  the  analysis  did  not  support  the  postulated  continuous-time 
external  feedback  control  mechanism.  The  time  internal  from  the  start  of  the 
saccade  to  the  time  at  peak  velocity  for  subjects  2  and  3  did  show  marked 
variation  ranging  from  21  msee  to  8  msec,  however.  Most  of  the  values  were 
clustered  about  evenly  into  two  grouping  16  msec  and  10  msec.  These  time 
interval  results  are  plotted  in  Fig.  5.  It  is  interesting  to  note  that  the  5° 
movements  elicited  under  natural  conditions  were  always  in  the  larger  time 
interval  grouping  of  16  msec  for  each  of  the  three  subjects.  It  is  speculated 
that  if  there  were  a  cont inuous-t ime  external  feedback  control  mechanism, 
then  both  the  time  interval  and  the  final  eye  displacement  would  have  been 
much  larger  than  under  normal  conditions. 

The  results  of  the  analysis  of  the  data  on  the  three  subjects  using  the 
system  identification  technique  confirmed  the  previous  findings. 
Additionally,  estimates  of  the  pulse  magnitude  of  the  agonist  active  state 
tension  F^p  wore  clustered  about  evenly  into  two  groupings,  corresponding  to 
the  previous  groupings,  0.95  NTS.  and  0.85  NTS.  These  results  are  plotted  in 
Fig.  6A.  The  estimated  pulse  magnitude  for  the  agonist  active  state  tension 
as  a  function  of  displacement  for  the  natural  saccades  are  illustrated  in  Fig. 


Fig.  (A.  The  estimated  ▼aloes  of  the  poise  magnitude  of  the  agonist 
aotWe  state  tension  ts  absolnte  ▼aloe  of  displacement  for 
aobjeot  3  dnring  the  stage  2  experiment. 


I 


daring  tha  ataga  1  axparinant 


Pulse  Magnitude 


6B.  One  pronounced  feature  evident  from  these  two  graphs  is  that  there  does 
not  appear  to  he  a  strong  relationship  between  pulse  magnitude  and 
displacement.  In  fact  the  variation  in  F^p  seen  in  Fig.  6A  is  approximately 
the  same  size  variation  seen  in  Fig.  6B.  Apparently  the  magnitude  of  the 
agonist  pulse  is  maximum  regardless  of  the  amplitude  of  the  saccade,  and  only 
the  duration  of  the  agonist  pulse  affects  the  size  of  the  saccade.  Under 
these  conditions,  the  eyeball  is  driven  to  its  destination  in  minimim  time  for 
saccades  of  all  sizes.  (See  Appendix  A3  for  a  more  complete  discussion  of 
this  time  optimal  model).  (21) 

VI.  KECOMXENTATIONS 


The  results  obtained  during  the  performance  of  this  project  indicate 
that  the  saccadic  eye  movement  system  operates  under  a  time  optimal  control 
strategy  without  continuous-time  external  feedback  sensory  information.  The 
implication  that  the  eyeball  is  directed  to  its  desination  in  minimum  time  for 
saccades  of  all  sizes  certainly  deserves  additional  attention  considering  the 
small  number  of  subjects  tested.  This  new  hypothesis  certainly  should  be 
investigated  further  in  light  of  the  current  theory  on  the  neuronal  control  of 
saccadic  eye  movements.  Particular  attention  should  be  paid  to  the  large 
variation  seen  with  saccades  of  the  same  size.  Furthermore.  the  importance 
of  the  second  phase  of  the  saccade,  that  is,  when  the  agonist  pulse  is  falling 
to  its  steady  state  value,  must  be  understood  in  relation  to  the  time  optimal 
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abstract 


The  objective  of  this  study  was  the  modification  of  a  linear 
homeomorphic  horizontal  saccadic  eye  movement  model  to  a  direct 
programming  state-space  representation  through  Laplace  variable 
analysis  about  the  operating  point  or  initial  eye  position.  The 
lateral  and  medial  rectus  muscle  of  each  eye  is  modeled  as  a 
parallel  combination  of  an  active  state  tension  generator  with  a 
viscosity  and  elastic  element,  connected  to  a  series  elastic 
element.  The  eyeball  is  modeled  as  a  sphere  connected  to  a 
viscosity  and  elastic  element.  Each  of  these  elements  is  assumed 


INTRODUCTION 


Models  of  oculomotor  function  are  important  in  the  develop¬ 
ment  of  clinically  useful  diagnostic  tools  and  in  understanding 
the  neurophysiology  of  eye  movements  [1]  -  [8).  The  complexity 
of  these  models  and  their  improved  correlation  with  physiological 
data  has  increased  since  Westheimer  first  presented  a  model  of 
saccadic  eye  movements  in  1954  [9]  .  Recently,  a  sixth-order 
linear  homeomorphic  horizontal  saccadic  eye  movement  model  was 
developed  that  provided  an  excellent  match  between  the  model 
predictions  and  the  data  [10],  However,  the  complexity  of  the 
model  has  hampered  its  widespread  use  in  clinical  studies  of 
oculomotor  function.  Investigators  have  utilized  zero-order  or 
second-order  oculomotor  plant  functions  with  reduced  precision 
and  little  or  no  theoretical  considerations  11]  -  [3]  [5]  [6] 
[11],  This  paper  presents  a  modification  of  the  linear  homeomor¬ 
phic  saccadic  eye  movement  model.  The  simplicity  of  the  derived 
model  makes  it  ideal  for  use  in  the  development  of  more  sensitive 
tests  of  oculomotor  pathology  and  in  the  description  of  normal 
oculomotor  function. 


MODEL  OP  OCULOMOTOR  PLANT 

Figure  1  illustrates  the  mechanical  components  of  the 
oculomotor  plant  for  horizontal  eye  movements,  the  lateral  and 
medial  rectus  muscle,  and  the  eyeball.  The  agonist  muscle  is 


modeled  as  a  parallel  combination  of  an  active  state  generator 
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Fag,  viscosity  element  BAG  and  elastic  element  K^,  connected  to 
a  series  elastic  element  Kgg.  The  antagonist  muscle  is  similarly 
modeled  as  a  parallel  combination  of  an  active  state  generator 
FANT'  viscosity  element  BANT  and  elastic  element  KLT,  connected 
to  a  series  elastic  element  KgE.  The  eyeball  is  modeled  as  a 
sphere  with  moment  of  inertia  Jp  connected  to  viscosity  element 
Bp  and  elastic  element  Kp.  The  passive  elasticity  for  each 
muscle  is  included  in  spring  Kp  for  ease  in  analysis.  Each  of 
the  elements  defined  in  the  model  of  the  oculomotor  plant  is 
ideal  and  linear. 

The  model  presented  in  Fig.  1  is  consistent  with  the  linear 
homeomorphic  model  except  for  the  treatment  of  initial  condi¬ 
tions.  The  muscles  are  shown  to  be  extended  from  equilibrium,  a 
position  of  rest,  at  the  primary  position  (looking  straight 
ahead)  consistent  with  physiological  evidence  [12] .  The  average 
length  of  the  rectus  muscle  at  the  primary  position  is 
approximately  40  mm,  and  at  the  equilibrium  position  is 
approximately  37  mm  [12].  0  is  the  angle  the  eyeball  is  deviated 
from  the  primary  position,  and  variable  x  is  the  length  of  the 
arc  traversed.  When  the  eye  is  at  the  primary  position,  both  0 
and  x  are  equal  to  zero.  Variables  through  x^  describe  the 
displacement  from  equilibrium  for  the  stiffness  elements  in  each 
muscle.  Values  Xp^  through  Xp^  are  the  displacements  from  equal- 
ibrium  for  each  of  the  variables  x^  through  x^  at  the  primary 
position.  The  total  extension  of  the  muscle  from  equilibrium  at 
the  primary  position  is  Xp^  plus  Xp2  or  Xp^  plus  Xp^  which 
equals  approximately  3mm  [12] .  It  is  assumed  that  the  lateral 
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and  medial  rectus  muscles  are  identical,  such  that  x^  equals 
Xp4  and  Xp3  equals  Xp2*  The  following  relationships  (evident  from 
Fig.  1)  exist  among  variables  x^,  x4,  ft,  and  x. 

xi  “  *  "  x  -,  (1) 


x4  =  X  +  xp4,  (2) 

and  «  =  57.296  x/r,  (3) 

where  r  equals  the  radius  of  the  eyeball  with  a  value  of  approx¬ 
imately  11  mm.  The  initial  condition  is  assigned  to  variables 


x4»  and  ®  through  x;  that  is,  variable  x  at  time  0,  x(0)  ,  and 
equations  1  through  3.  Initial  conditions  are  also  assigned  to 
variables  x2  and  x3  as  x2(0)  and  x3^°^»  and  to  the  agonist  and 
antagonist  active  state  tensions  as  FAG(®)  and  PANT^  * 

A  set  of  three  equations  is  written  to  describe  the  oculo¬ 
motor  plant  by  summing  the  forces  acting  at  junctions  2  and  3, 
the  equilibrium  positions  for  x2  and  x^,  and  the  torques  acting 
on  the  eyeball. 


FAG  “  KLTx2  + 

BAG*2 

+  KSE^X2  “  xl^ ' 

(4) 

and 

FANT  =  KSEtx4 

-  x3> 

KLTX3  ~  BANTX3  * 

•  •  • 

(5) 

rKSE(x2  +  x3  -  x!  - 

X 

II 

J  ft  +  B  ft  +  K  0. 

P  P  P 

(6) 

Equations  4  and  5  differ  from  those  of  Bahill  ££.  al.  with  regard 
to  the  viscosity  terms,  a  constant  KSE/(KLT  +  KgE)  is  missing  in 
the  viscosity  term  in  their  equations  14  and  15  [10J.  Equations 
4  through  6  could  be  simplified  and  written  as  one  equation  by 
taking  the  Laplace  transform  of  each  equation  and  eliminating  all 
variables  except  x.  However,  the  initial  conditions  must  be 
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and 


K  =  57.296  Kp/r  . 

Note  that  variables  x^  and  x4  are  eliminated  from  the  above  equa- 

A 

tions  since  they  are  both  equal  to  x,  and  variable  0  through 
equation  3. 

Equations  7  through  9  are  written  as  one  equation  by  taking 
the  Laplace  transform  of  each  equation  and  eliminating  variables 
x2  and  x3  yielding 

kse*sbant  +  kst*l*fag*  ”  kse*sbag  +  kst)l{fant} 

(C4s4  +  C3S3  +  C2S2  +  CXS  +  C0)  X(S) ,  (10) 

where  C4  = 

C3  =  jkst(bag  +  bant*  +  bbantbag/ 

C2  =  jkst  +  bkst*bag  +  bant*  +  bagbant^k  +  2kse* ' 

C1  =  bkst  +  *bag  +  bant*  <kkst  +  2ksekst  ~  kse*  ' 
and  C0  =  KKst  +  2kseKSTKLT* 

Transforming  equation  10  back  into  the  time  domain  yields 

BA3FANT^  =  C4*  +  ^3*  +  ^2*  +  *1*  +  S)** 

This  equation  describes  the  trajectory  of  the  eyeball  during  a 
horizontal  saccadic  movement  from  the  initial  eye  position.  By 
introducing  the  initial  eye  position  and  active  state  tension, 
and  noting  that 


kse*kst(fag  fant)  +  bantfag 
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P  c  P 

ANT  ANT' 


we  have 


*WKSTtPAG  ~  FANT^  +  BANTPAG“BAGPANT)  “  C4X  +  C3X*  +  ^2*  +  °1*  +  C0X*  *12) 

Since  data  are  available  for  the  angle  of  the  eyeball  with 
respect  to  the  primary  position,  the  angle  0  is  used  instead  of  x 
which  yields 


6lKSTtFAG“FANT^+BANTFAG“BAGFANT^SS  9  +  P39  +  P29  +  Pl®  +  P09/  (13) 


where 


6  =  57.296KSE/rC4, 

P3  =  C3^CA 

P2  =  C2^C4' 

P1  *  Cj/C4' 

and  Pq  =  CQ/C4 . 

State  equations  are  usually  written  for  ease  in  simulating 
saccadic  eye  movements  with  programs  available  for  the  digital 
computer  such  as  the  Continuous  System  Modeling  Program.  One  set 
of  state  equations  is  given  below  in  terms  of  state  variables  8^ 
through  ©4,  where 

©2  =  8  =  angular  position, 

02=02=  angular  velocity, 

©2  =  0£  =  angular  acceleration, 
and  ©4  =  83  =  angular  jerk, 

with  state  equations 

=  e2, 

9,  -  e„ 
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*3  "  V 


04  =  "P0®1  “  Pl02“P2e3~P364+  (KSTtFAG“FANT*+BANTFAG“BAGFANT* * 

The  agonist  and  antagonist  active  state  tension  terms  are  not 
included  as  state  variables  since  they  do  not  describe  the  energy 
state  of  the  system. 


DISCUSSION 

The  objective  of  this  study  was  to  modify  the  linear  homeo- 
morphic  model  to  better  understand  the  complex  mechanisms  under¬ 
lying  saccadic  eye  movements,  and  to  derive  a  simpler  repre¬ 
sentation  of  the  model  in  an  effort  to  increase  its  potential 
value  in  evaluating  normal  and  pathological  oculomotor  function. 
We  feel  we  have  accomplished  this  objective  as  illustrated  in  the 


previous  section. 


Since  no  additional  assumptions  are 


introduced,  the  excellent  match  between  model  prediction  and  the 
data  reported  by  Bahill  jei  al,  remains  intact  110] .  The  elimina¬ 
tion  of  variables  x2  and  x^  from  the  linear  homeomorphic  model 
leads  to  a  simpler  direct  programming  state-space  representation. 
Note  that  the  state  variables,  velocity,  acceleration  and  jerk, 
are  estimable  from  the  measured  position  of  the  eyeball  either 
from  analog  electronic  processing  or  through  computational 
techniques  such  as  the  two  point  central  difference  method  [13] 
114].  Human  physiological  data  is  available  for  estimating  each 
parameter  in  the  linear  homeomorphic  model  [12]  [15].  While  the 


tension  produced  by  each  muscle  is  experimentally  known  115}, 


the  actual  active  state  tensions  (muscle  forces),  P^G  and  FANT, 
and  the  initial  conditions  on  variables  through  x^  are 
unknown.  Since  the  initial  conditions  for  variables  x^  through 
x^  are  not  needed  in  the  modified  linear  homeomorphic  model,  only 
the  saccadic  muscle  forces  need  be  specified.  One  typically 
assumes  an  input  waveform  for  the  agonist  and  antagonist  active 
state  tensions  that  is  a  simple  function  of  the  observed 
motoneuronal  signal  116] [17} {18} . 

Presented  in  Fig.  2  is  a  block  diagram  of  the  saccadic  eye 
movement  system  which  functions  as  one  element  within  the  overall 
system  controlling  all  oculomotor  movements.  The  saccadic  eye 
movement  system  consists  of  three  components;  the  controller, 
that  is,  the  agonist  and  antagonist  active  state  tension,  the 
oculomotor  plant,  and  the  feedback  element  H.  Since  no  assump¬ 
tions  concerning  the  saccadic  active  state  tensions  are  required 
in  the  derivation  of  the  modified  linear  homeomorphic  model,  any 
type  of  control  signal  may  be  utilized.  For  example,  investi¬ 
gators  may  utilize  this  model  when  examining  alternatives  to  the 
pulse-step  muscle  force.  The  oculomotor  plant,  indicated  by  the 
section  within  the  dashed  lines  in  Fig.  2,  is  based  on  equation 
10.  The  four  real  eigenvalues  as  determined  from  the 
characteristic  equation  for  the  parameter  values  given  by  Bahill 
et  al .  are  -17,  -75,  -171,  and  -1307  [10]. 

The  mechanical  components  of  the  model  are  connected  such 
that  the  rate  of  change  of  the  muscle  forces  together  with  the 
actual  muscle  forces  drive  the  eyeball  to  its  final  position. 
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Consider  the  low-pass  filtered  pulse-step  active  state  tension 
signals  presented  in  Fig.  3.  During  the  initial  pulse  phase  of 
the  trajectory,  the  agonist  active  state  tension  drives  the 
eyeball  and  the  antagonist  active  state  tension  restrains  the 
eyeball.  Since  the  antagonist  active  state  tension  falls  to  zero 
during  this  interval,  then  the  rate  of  change  antagonist  active 
state  tension,  is  negative  and  as  indicated  by  the  equation 
13,  acts  to  drive  the  eyeball  along  with  the  agonist  force,  F*-, 
and  the  agonist  rate  of  change  force,  F^G .  As  a  result,  a 
greater  force  propels  the  eyeball  than  if  only  the  agonist  and 
antagonist  active  state  tensions  are  used  as  shown  in  the  plant 
input  curves  of  Fig.  3.  This  result  is  in  agreement  with  the 
tension  data  recorded  for  horizontal  saccadic  eye  movements 
[15JI19],  After  the  pulse  phase  of  the  trajectory, ‘the  agonist 
active  state  tension  falls  to  a  steady  state  value  which  results 
in  a  negative  rate  of  change  agonist  active  state  tension,  FAG* 

Thus,  according  to  equation  13,  the  agonist  rate  of  change  force, 

0 

FAG'  acts  to  restrain  the  eyeball  with  the  antagonist  force, 
FANTf  and  the  anta90nist  rate  of  change  force,  FANT*  As  shown 
for  the  plant  input  curves  in  Fig.  3,  the  resultant  forces  act 
like  a  dynamic  brake  which  is  reduced  to  zero  as  the  eyeball 
reaches  its  destination. 

In  summary,  a  modification  of  a  linear  homeomorphic  model 
for  horizontal  saccadic  eye  movements  is  presented  which  allows 
one  to  estimate  or  directly  measure  relevant  parameters  and 
initial  conditions  from  physiological  data.  Since  no  assumptions 
concerning  the  saccadic  active  state  tensions  are  used  in  the 


10 

derivation  of  the  model,  any  type  of  control  input  can  be  used  in 
evaluating  the  saccadic  system. 
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The  mechanical  components  for  the  oculomotor  plant 
at  the  primary  position. 


A  block  diagram  of  the  oculomotor  plant.  The 
section  of  the  diagram  within  the  dashed  lines  is 
the  modified  linear  homeomorphic  model.  The 
feedback  element  H  is  unity. 


Transformation  of  the  active  state  tensions  by  the 
numerator  terms  of  the  oculomtor  plant.  One 
diagram  shows  the  assumed  low-pass  filtered 
pulse-step  waveform  of  the  agonist  and  antagonist 
active  state  tensions  initiated  at  100  msec.  The 
antagonist  pulse  circumscribes  the  agonist  pulse,  by 
3  msec  on  each  end.  Parameter  values  are  from 
Bahill  fit  al.  for  simulating  a  10  degree  horizontal 
saccadic  eye  movement  [101.  The  other  diagram 
shows  a  partial  and  full  transformation  of  the 
active  state  tensions  by  the  numerator  terms  of  the 
oculomotor  plant. 


Two-point  central  difference  analysis  results  on  three  subjects 
for  natnral  and  modified  target  trajectory  saccade  data. 

List  of  Symbols 

Saeeade  -  distance  the  target  moved 

Left  and  Sight  -  the  distance  the  left  and  right  eye  moved  in 
DisplaeeMnt  response  to  the  target  displacement  (degrees) 

I.C.  -  initial  position  of  the  eyeball  (degrees) 

TD  -  time  from  the  start  of  the  saccade  to  the  time  at  peak 
velocity  (msec) 

DUE  -  duration  of  the  saccade  (msec) 

TP  -  time  at  which  the  saccade  started  (msec) 

P?  -  peak  velocity  of  the  saccadic  eye  movement  (degrees/second) 
DIST  -  eye  position  at  the  time  of  peak  velocity  (degrees) 


Normal  Saccade  Data  for  Subject  1 


Normal  Saccade  Data  for  Subject  3 


OPTIMAL  CONTROL  OP  SACCADIC  EYE  MOVEMENTS 

John  D.  Enderle 
Division  of  Bioengineering 

Department  of  Electrical  and  Electronics  Engineering 
North  Dakota  State  University 
Fargo,  North  Dakota  58105 
and 

James  W.  Wolfe  and  J.  Terry  Yates 
Neurosciences  Function 
USAF  School  of  Aerospace  Medicine 
Brooks  Air  Force  Base 
San  Antonio,  Texas  78235 


Proceedings  of  the  Twenty-First  Annual  Rocky  Mountain 
Bioengineering  Symposium,  Boulder,  Colorado,  April  1984 


ABSTRACT 

A  new  theory  describing  the  optimal  control  of  saccadic 
eye  movements  is  proposed  based  on  physiological 
considerations.  The  lateral  and  medial  rectus  muscle  of 
each  eye  is  modeled  as  a  Voigt  element  connected  to  a  series 
elastic  element.  The  eyeball  is  modeled  as  a  sphere  con¬ 
nected  to  a  viscosity  and  elastic  element.  Each  of  these 
elements  is  assumed  to  be  ideal  and  linear.  The  active 
state  tension  for  each  muscle  is  modeled  by  a  low-pass 
filtered  pulse-step  waveform.  The  magnitude  of  the  agonist 
pulse  is  assumed  to  be  a  maximum  constant  for  saccades  of 
all  sizes  and  only  the  duration  of  the  agonist  pulse 
affects  the  size  of  the  saccade.  The  antagonist  muscle  is 
assumed  to  be  completely  inhibited  during  the  period  of 
maximum  stimulation  for  the  agonist  muscle.  Under  these 
assumptions,  the  eyeball  is  directed  to  its  destination  in 
minimum  time  for  saccades  of  all  sizes.  Horizontal  saccadic 
eye  movements  were  recorded  from  infrared  reflection  signals 
off  the  anterior  surface  of  the  cornea  and  then  digitized. 
Parameter  estimates  are  calculated  for  the  model  by  using  a 
conjugate  gradient  search  program  which  minimizes  the 
integral  of  the  absolute  value  of  the  error  squared  between 
the  model  and  the  data.  The  predictions  of  the  model  are 
shown  to  be  in  good  agreement  with  the  data. 


INTRODUCTION 


Saccadic  eye  movements  are  among  the  fastest  voluntary 
muscle  movements  the  human  body  is  capable  of  producing. 
Although  the  purpose  for  such  eye  movement  is  obvious,  that 
is,  to  quickly  redirect  the  eyeball  on  the  target,  the 
neuronal  control  strategy  is  not.  For  instance,  does  the 
word  "quickly"  in  the  previous  sentence  imply  the  most  rapid 
movement  possible,  or  simply  a  fast  as  opposed  to  a  slow 
movement?  One  group  of  researchers  investigating  the 
neuronal  control  strategy  concluded  that  each  saccadic  eye 
movement  is  optimally  driven  by  a  low-pass  filtered 
pulse-step  muscle  force  (active  state  tension)  to  achieve 
the  final  eye  position  in  minimum  time  (1) (2) .  Their 
results,  however,  failed  when  applied  as  a  consistent 
neuronal  control  mechanism  for  all  saccadic  eye  movements 
independent  of  the  magnitude  of  the  saccade.  In  order  to 
reach  a  destination  in  minimum  time,  the  system  should  be 
fully  stimulated,  and  then  fully  braked,  regardless  of  the 
distance  traversed.  Under  this  control  strategy,  the  only 
variable  affecting  saccade  magnitude  is  the  length  of  time 
that  the  system  is  fully  stimulated.  Since  they  assumed 


that  saccade  magnitude  is  both  a  function  of  the  amplitude 
and  duration  of  the  saccadic  innervation  signal,  and  not  of 


duration  alone/  it  is  not  an  obvious  minimum  time 
controller . 

In  order  to  detail  the  neuronal  control  strategy,  it  is 
necessary  to  understand  the  effect  of  the  measurable 
saccadic  innervation  signal  on  the  oculomotor  plant  and  the 
resultant  eyeball  response.  Collins  states  that  the 
magnitude  of  each  saccade  is  determined  by  the  amplitude  and 
duration  of  the  saccadic  innervation  signal  (3). 
Specifically,  he  determined  a  logarithmic  innervation 
amplitude  to  saccade  magnitude  relationship.  Zee  and  his 
co-workers  assumed  that  the  amplitude  and  duration  of  the 
saccadic  innervation  signal  is  automatically  controlled  by  a 
local  feedback  loop  (4).  The  saccadic  innervation  signal  is 
driven  by  the  difference  between  the  internal  representation 
of  the  present  eye  position  and  the  desired  eye  position. 
Note  that  these  authors  still  hypothesize  a  pulse-step 
innervation  signal,  but  the  pulse  size  and  duration  are 
based  on  a  nonlinear  velocity  function.  If  these  authors 
are  correct  in  their  assertion,  that  is,  saccade  magnitude 
is  a  function  of  both  the  duration  and  the  amplitude  of  the 
saccadic  innervation  signal  (and  not  duration  alone) ,  then 
the  neuronal  control  does  not  operate  with  an  apparent 
minimum  time  strategy. 

Recently,  Robinson  presented  data  that  seems  to 
contradict  the  previous  saccadic  innervation  relationship  of 


amplitude  and  duration  with  saccade  magnitude  (5). 
Innervation  seems  to  peak  at  the  same  height  and  maintain  a 
constant  but  lower  level  throughout  the  saccade  for  the 
agonist  muscle  for  any  size  retinal  error  (see  Fig.  4  in 
Robinson) .  Apparently,  only  the  duration  of  the  agonist 
pulse  is  a  function  of  the  target  displacement.  A  possible 
explanation  for  this  disparity  with  other  investigators 
recordings  may  lie  in  the  type  of  motoneuron  from  which  the 
recordings  were  made  (3).  If,  in  fact,  Robinson’s  saccadic 
innervation  signals  are  a  correct  representation,  then  a 
neuronal  control  strategy  of  minimum  time  still  does  not 
seem  appropriate  because  of  the  higher  initial  firing  rates. 

While  investigators  have  been  able  to  record  the 
innervation  signal  from  several  types  of  motoneurons  that 
drive  the  eyeball  during  a  saccade,  they  have  not  been  able 
to  directly  measure  the  muscle  force  responsible  for  this 
movement.  Collins  and  his  co-workers  have  measured  the 
tension,  in  vivo,  at  the  muscle  .tendon  during  unrestrained 
human  eye  movement  using  a  miniature  "C"  gauge  force  trans¬ 
ducer  (6).  The  muscle  force,  however,  is  distributed 
throughout  the  muscle  and  cannot  be  directly  measured  since 
it  is  modified  by  the  visoelasticity  of  the  muscle.  Very 
little  is  known  about  the  dynamic  muscle  forces  generated  in 
the  antagonist-agonist  pair  during  a  saccade  and  their 
relationship  to  the  saccadic  innervation  signal.  Naturally, 


under  static  conditions  during  fixation,  the  muscle  forces 
are  proportional  to  the  innervation  signal,  the  constant 
firing  frequency  produces  a  constant  force.  During  a 
saccade,  the  agonist  muscle  force  changes  rapidly,  rising  to 
a  new  level  approximately  ten-fold  higher  than  fixation  in 
milliseconds,  and  then  falling  to  the  new  fixation  level. 
The  proportional  relationship  that  exists  for  innervation 
and  muscle  force  during  fixation  is  not  valid  during  a 
saccade  due  to  the  effect  of  saturation  and  filtering  on  the 
innervation  signal.  Thus,  the  exact  shape  of  the  muscle 
waveform  is  uncertain.  It  has  been  observed  that  the  muscle 
force  lags  behind  the  saccadic  innervation  signal  by  8  msec 
due  to  diffusion  and  resequester.ing  of  Ca++  ions  in  the 
muscle  and  is  typically  modeled  by  low-pass  filtering  of  the 
innervation  signal  (5).  While  very  little  is  known  about 
the  activation  and  deactivation  time  constant  due  to  lack  of 
in  vivo  testing,  Bahill  has  estimated  their  values  to  be 
between  .2  and  13  msec  based  on  the  rise  of  the  isometric 
force  during  electrical  stimulation  (7).  When  muscles  are 
artificially  stimulated,  they  do  not  develop  additional 
tension  when  the  stimulus  frequency  is  above  200  Hz  (5). 
The  normal  firing  frequency  for  the  agonist  muscle  averages 
approximately  800  Hz  during  the  saccade  (5).  One  explana¬ 
tion  for  why  motoneurons  fire  at  rates  beyond  which  the 
muscle  can  respond  is  that  both  the  antagonist  and  agonist 


rate  of  change  force  significantly  contributes  to  driving 
the  eyeball  to  its  destination  (5) (8) (9) .  Since  saccades  of 
all  sizes  fire  well  above  200  Hz  for  the  initial  pulse  phase 
of  the  trajectory,  then  are  the  saccades  all  driven  by  the 
same  magnitude  agonist  muscle  force  pulse?  Furthermore,  if 
the  magnitude  of  the  agonist  pulse  is  the  same  for  any  size 
saccade,  then  the  only  variable  affecting  the  saccadic  size 
is  the  duration  of  the  agonist  pulse.  Under  these 
conditions,  the  eyeball  is  driven  to  its  destination  in 
minimum  time  for  saccades  of  all  sizes. 

In  this  paper  we  present  a  quantitative  analysis  of 
saccadic  eye  movements  to  support  the  hypothesis  that  the 
saccadic  neuronal  control  mechanism  operates  to  achieve 
final  eye  position  in  minimum  time.  The  concepts  underlying 
this  hypothesis  are:  (1)  the  active  state  tension  for  each 
muscle  is  modeled  by  a  low-pass  filtered  pulse-step  waveform 
in  which  the  magnitude  of  the  agonist  pulse  is  a  maximum 
constant  for  saccades  of  all  sizes,  and  (2)  only  the 
duration  of  the  agonist  pulse  affects  the  size  of  the 
saccade.  The  antagonist  muscle  is  assumed  to  be  completely 
inhibited  during  the  period  of  maximum  stimulation  for  the 
agonist  muscle. 


Data  were  collected  after  seating  a  subject  before  a 
target  display  of  nine  small  red  light  emitting  diodes 
(LED),  each  separated  by  five  degrees.  The  subject's  head 
was  restrained  by  a  bite-bar.  The  subject  was  instructed  to 
follow  the  "jumping"  target  which  moved  from  the  center 
position  to  any  of  the  other  LED's  and  then  returned  to  the 
center  position.  The  order  of  the  target  positions  as  well 
as  the  time  interval  between  displacements  was  randomized. 
Data  were  only  recorded  for  the  initial  displacement  from 
the  center  position. 

Horizontal  eye  movements  were  recorded  from  each  eye 
using  an  infrared  reflection  signal  from  the  anterior 
surface  of  the  cornia-scleral  interface  with  instrumentation 
described  by  Engelken  et  al.  (10)  .  Signals  for  both  eyes 
tracking  were  digitized  using  the  analog/digitial  converter 
of  a  DBCLAB  PDP  11/34  computer  and  stored  in  disk  memory. 
These  signals  were  sampled  at  a  rate  of  1000  samples  per 
second  for  one-half  second. 

MODEL  OF  THE  OCULOMOTOR  PLANT 

The  diagram  in  Fig.  1  illustrates  the  mechanical 
components  of  the  oculomotor  plant  for  horizontal  eye 
movements.  The  lateral  and  medial  rectus  muscle  of  each 
eye  is  modeled  as  a  Voigt  element  connected  to  a  series 


elastic  element  KSE*  A  Voigt  element  for  the  agonist  muscle 
consists  of  an  ideal  active  state  tension  generator  PA^  in 
parallel  with  viscolsity  element  and  length  tension 

spring  KLT.  The  eyeball  is  modeled  as  sphere  with  moment  of 
inertia  Jp  connected  to  viscosity  element  Bp  and  elastic 
element  Kp.  The  passive  elasticity  for  each  muscle  is 
included  in  spring  Kp  for  ease  in  analysis.  Each  of  the 
elements  defined  in  the  model  of  the  oculomotor  plant  is 
ideal  and  linear.  0  is  the  angle  the  eyeball  is  deviated 
from  the  primary  position  (looking  straight  ahead),  and 
variable  x  is  the  length  of  the  arc  traversed.  When  the  eye 
is  at  the  primary  position,  both  0  and  x  are  equal  to  zero. 
Variables  x^  through  x^  describe  the  displacement  from 
equalibrium  for  the  stiffness  elements  in  each  muscle. 
Values  xpl  through  xp4  are  the  displacements  from 

equalibrium  for  each  of  the  variables  x^  through  x4  at  the 
primary  position.  The  total  extension  of  the  muscle  from 
equilibrium  at  the  primary  position  is  Xp^  plus  Xp2  or 
Xp^  plus  Xp4  which  equals  approximately  3mm  (11)  .  It  is 
assumed  that  the  lateral  and  medial  rectus  muscles  are 
identical,  such  that  xp^  following  fourth-order  model  was 
derived  for  the  oculomotor  system  shown  in  Fig.  1  by 
Enderle,  V7olfe,  and  Yates  (9). 


{kst (fag  fant^  +bantfag  bagfakt* 
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RESULTS 

Representative  saccadic  eye  movement  data  and  model 
simulation  results  are  shown  in  Fig.  2.  A  similar  close 
agreement  between  the  predicted  saccadic  eye  response  and 
the  data  is  seen  for  input  target  displacements  from  -20  to 
23  degrees.  The  simulation  results  are  calculated  using  the 
IBM  Continuous  System  Modeling  Program  for  the  oculomotor 
model  in  equation  1,  and  the  following  low-pass  filtered 
pulse-step  signals  for  the  agonist  and  antagonist  active 
state  tensions. 

=  FAGOU  ^ ^  +  ^FAGP*”FAGO^  /Tac^  ^  U 

^FAGP~FAGS^  (1"exP{~{t“tp-td)/Tde^u^t-tp“td^ ' 


FANT(t,=  FANTO  PANTO  ^1-exPt  ~  (t-tp)  /  T  de})u(t-tp)  + 

FANTS (1“exp^" {t_tp"tdi  Ade^J u tt_tp"td)  ' 


where 

FAGO 

fAGP 

fags 

fanto 

fants 


ac 

and 


initial  magnitude  of  the  agonist  active  state 
tension 

pulse  magnitude  of  the  agonist  active  state 
tension, 

step  magnitude  of  the  agonist  active  state 
tension, 

initial  magnitude  of  the  antagonist  active  state 
tension, 

step  magnitude  of  the  antagonist  active  state 
tension, 

latent  period,  the  time  interval  between  the  target 
movement  and  the  start  of  the  eye  movement 
duration  of  the  agonist  pulse  active  state  tension 
activation  time  constant, 


r  de  =  deactivation  time  constant. 


Quantitative  analysis  of  the  data  was  conducted  in  two 
stages.  During  the  first  stage,  all  model  and  force 
parameters  are  estimated  with  a  conjugate  gradient  search 
program  which  minimizes  the  integral  of  the  absolute  value 


stages.  During  the  first  stage,  all  model  and  force 
parameters  are  estimated  with  a  conjugate  gradient  search 
program  which  minimizes  the  integral  of  the  absolute  value 
of  the  error  squared  between  the  frequency  response  and  the 
data  for  the  first  saccadic  eye  movement  (12) . 
Physiological  data  is  used  to  estimate  all  of  the  initial 
parameter  estimates  (6) (11) .  During  the  second  stage,  all 
subsequent  saccadic  eye  movements  are  analyzed  using  the 
conjugate  gradient  search  program  in  which  only  the  latent 
period  and  the  duration  of  the  agonist  pulse  are  optimally 
adjusted.  All  other  model  and  force  parameter  values  are 
set  equal  to  those  values  computed  during  the  first  stage. 

DISCUSSION 

The  objective  of  this  study  was  to  investigate  the 
hypothesis  that  saccadic  eye  movements  are  generated  under  a 
time  optimal  neuronal  control  strategy  for  saccades  of  all 
sizes.  Specifically,  we  postulate  that  the  active  state 
tension  for  each  muscle  is  modeled  by  a  low-pass  filtered 
pulse-step  waveform  in  which  the  magnitude  of  the  agonist 
pulse  is  a  maximum  constant  for  saccades  of  all  sizes  and 
that  only  the  duration  of  the  agonist  pulse  affects  the  size 
of  the  saccade.  The  antagonist  muscle  is  assumed  to  be 
completely  inhibited  during  the  period  of  maximum 
stimulation  for  the  agonist  muscle.  Thus,  the  eyeball  is 


driven  to  its  final  position  in  minimum  time  for  saccades  of 
all  sizes.  The  time  optimal  hypothesis  is  supported  by  the 
simulation  results  which  closely  match  the  measured  saccadic 
eye  movement  data  for  target  displacements  from  -20  to  20 
degrees.  A  time  optimal  neuronal  control  strategy  is  also 
supported  through  physiological  evidence.  Experiments 
involving  artificial  stimulation  of  oculomotor  muscles  above 
200  Hz  indicate  that  such  muscles  do  not  develop  additional 
tension  (5) .  Since  the  agonist  firing  rate  for  saccadic  eye 
movements  is  much  above  200  Hz,  then  the  same  magnitude 
agonist  pulse  must  drive  the  eyeball  regardless  of  saccade 
size.  Experimental  results  indicate  a  clear  relationship 
between  the  duration  of  the  agonist  saccadic  innervation 
signal  (and  corresponding  active  state  tension)  and  the  size 
of  the  saccade. 

It  is  well  known  in  optimal  control  theory  that  in 
order  to  transfer  a  system  from  point  A  to  point  B  in 
minimum  time,  the  system  should  be  maximally  accelerated  for 
one-half  the  time  interval,  and  then  maximally  decelerated 
for  \the  remainder  of  the  time  interval  (13)  .  Due  to  the 

v 

\ 

low-  pass  filtering  on  the  active  state  tension  signal  and 
the  unequal  activation  and  deactivation  time  constants,  the 
periods  of  acceleration  and  deceleration  are  unequal  for 
saccadic  eye  movements.  Evidence  indicates  that  during  the 
period  of  acceleration,  the  agonist  pulse  takes  on  a  maximum 


value  regardless  of  the  size  of  the  saccade.  During  the 
period  of  deceleration,  the  oculomotor  system  is  dynamically 
braked  by  the  antagonist  step  waveform  and  the  rate  of 
change  forces,  d  F^/dt  and  d  F^^/dt  (9).  Therefore,  the 
saccadic  eye  movement  system  is  in  fact,  a  time  optimal 
system. 
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LEFT  EYE  RESPONSE  IN  DEGREES 


Figure  2.  A  10°  Saccadic  Eye  Movement. 
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ABSTRACT 


A  simple  second-order  model  of  the  horizontal  saccadic  eye  move¬ 
ment  system  is  developed  from  physiological  considerations.  The 
lateral  and  medial  rectus  muscles  of  each  eye  are  modeled  as  a 
Voigt  element  connected  to  a  series  elastic  element.  The  eye¬ 
ball  is  modeled  as  a  sphere  connected  to  a  viscosity  element  and 
an  elastic  element.  Each  of  these  elements  are  assumed  to  be 
ideal  and  linear.  The  model  is  derived  so  that  the  initial 
conditions  for  the  elastic  elements  are  included  in  the  initial 
agonist  and  antagonist  muscle  forces  and  eye  position.  As  a 
result,  all  parameters  and  initial  conditions  can  be  estimated  or 
directly  measured  from  physiological  data.  The  derivation  is 
sufficiently  general  that  no  requirements  are  assumed  for  the 
saccadic  muscle  forces.  Close  agreement  between  the  predicted 
saccadic  eye  response  and  the  data  was  obtained  for  target  dis¬ 
placements  from  -20  to  +20  degrees  with  an  assumed  low-pass 
filtered  pulse-step  muscle  force. 


1.  INTRODUCTION 


Models  of  oculomotor  function  are  important  in  the  development  of 
clinically  useful  diagnostic  tools,  and  in  understanding  how  the 
human  brain  controls  eye  movement.  (l)-(8)  The  complexity  of 
oculomotor  models  and  their  correlation  with  physiological  evi¬ 
dence  has  increased  since  Westheimer  first  presented  a  model  of 
saccadic  eye  movements  in  1954.  (9)  Unfortunately,  the  complex¬ 
ity  of  these  models  have  precluded  their  widespread  use  in  most 
clinical  studies  of  oculomotor  function.  Investigators  have 
utilized  zero-order  or  second  order  plant  functions  with  reduced 
precision  and  little  or  no  theoretical  considerations.  (l)-(3) 
(5)  (6)  (10) 

This  paper  presents  a  linear  second-order  oculomotor  plant  model 
that  is  developed  from  physiological  data.  This  model  is  a 
modification  of  the  sixth-order  linear  homeomorphic  model  pre¬ 
sented  by  Bahill,  Latimer  and  Troost  which'  retains  the  excellent 
match  between  the  model  predictions  and  the  data.  (11)  The  sim¬ 
plicity  of  the  model  and  its  correlation  with  physiological  data 
make  it  ideal  for  use  in  the  development  of  more  sensitive  tests 
of  oculomotor  pathology  and  in  the  description  of  normal  oculo¬ 
motor  function. 

2.  MODEL  OP  OCULOMOTOR  PLANT 


Figure  1  illustrates  the  mechanical  components  of  the  oculomotor 
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Figure  1.  Mechanical  components  for  the  oculomotor  plant  at  the 
primary  position. 


plant  for  horizontal  saccadic  eye  movements.  The  lateral  and 
medial  rectus  muscles  of  each  eye  are  modeled  as  a  Voigt  element 
connected  to  a  series  elastic  element  KgE.  A  Voigt  element  for 
the  agonist  muscle  consists  of  an  ideal  force  generator  FftG  in 
parallel  with  viscosity  element  BAG  and  length  tension  spring 
K_  *.  The  eyeball  is  modeled  as  a  sphere  with  moment  of  inertia 
J  connected  to  viscosity  element  IS  and  elastic  element  K  .  The 


passive  elasticity  for  each  muscle  is  included  in  spring  for 
ease  in  analysis  1 

The  model  presented  in  Figure  1  is  consistent  with  the  linear 
homeomorphic  model  except  for  the  treatment  of  initial  condi- 
tions.  The  muscles  are  shown  to  be  extended  from  equilibrium  at 
the  primary  position  (looking  straight  ahead)  consistent  with 
physiologic  evidence.  (12)  Variables  Xj  through  x^  described  the 
displacement  from  equilibrium  for  the  stiffness  elements  in  each 
muscle.  Values  x^  through  Xp^  are  the  displacements  from  equil¬ 
ibrium  for  each  of  the  variables  x^  through  x^  at  the  primary 
position.  The  total  extension  of  the  muscle  from  equilibrium  at 
the  primary  position  is  xpl  plus  Xp2  or  xp3  plus  xp4  which  equals 
approximately  3mm.  (12)  0  is  the  angle  the  eyeball  is  rotated 
from  the  primary  position,  and  variable  x  is  the  length  of  the 
arc  traversed.  Each  of  the  elements  defined  in  Figure  1  is  ideal 
and  linear. 


A  set  of  three  linear  equations  is  written  to  describe  the 
oculomotor  plant  by  summing  the  forces  acting  at  junctions  2  and 

i 

3  and  the  torques  acting  on  the  eyeball. 


FAG  “  KLTX2  +  BAGX2  +  KSE(X2  ~  Xl) ' 
FANT  “  KSE(x4  X3J  KLTX3  "  BANTx3' 


(1) 

(2) 


and 

rKSB<x2  +  x3  '  X1  "  x4)  "  V  +  V  +  V'  (3> 

where  r  equals  the  radius  of  the  eyeball  with  a  value  of  approxi¬ 
mately  11  mm.  Equations  1  through  3,  which  described  the  trajec- 


tory  of  the  eyeball  during  horizontal  saccadic  movement,  are 
reduced  to  the  following  equation 


6<KST{FAG  “  PANT*  +  BANTPAG  “  BAGFANT)  * 

*0  +  P3*0  +  P  2§  +  PjO  +  Pq0  .  (4) 

Details  of  the  derivation  and  the  parameter  description  are 
contained  in  the  Appendix.  This  fourth-order  model  is  derived 
with  no  constraints  placed  on  the  type  of  force  used  to  drive  the 
eyeball  to  its  final  position,  such  as  a  pulse-step  force.  (13) 
An  assumption  has  been  made  that  there  is  a  constant  force  which 
acts  on  the  eyeball  until  it  moves.  This  notion  is  supported  by 
studies  of  the  discharge  patterns  of  motoneurons.  (14) 

It  is  possible  to  reduce  the  order  of  the  system  model  in  equa¬ 
tion  4  by  approximating  it  with  a  second-order  model.  (15) -(17) 
This  method  involves  minimizing  the  frequency  response  differ¬ 
ences  between  the  second-order  system  model  and  the  original 
system  model.  It  is  desired  to  maintain  the  numerator  of  the 
transfer  function  of  the  original  model  in  the  second-order  model 
for  ease  in  analysis.  This  is  accomplished  by  allowing  only  the 
denominator  of  the  second-order  model  to  be  subject  to  approxima¬ 
tion.  The  second-order  model  describing  horizontal  saccadic 
movement  is 

6  <Rst<fag  ~  FANT*  +  bantpag  ~  bagpant* 
d20  +  dj©  +  PQ0 ,  (5) 

where 

a2  -  (2p0  -  zpjPj  +  p£>1/2. 


(6) 


and 


dl  88  (“2P2P0  +  pi  +  2Pod2>1/2*  <7) 

The  eigenvalues  for  the  second-order  model  could  have  been 

estimated  using  the  two  dominant  poles  of  the  fourth-order  model. 

However,  this  would  have  resulted  in  less  precision  since  the 

« 

third  pole  affects  the  frequency  response  commencing  at  approxi¬ 
mately  30  radians  per  sec. 


3.  RESULTS 


Data  were  collected  after  seating  subject  before  a  target  display 
of  nine  small  red  light-emitting  diodes  (LED) ,  each  separated  by 
five  degrees.  The  subject's  head  was  restrained  by  a  bite-bar. 
The  subject  was  instructed  to  follow  the  "jumping"  target  which 
moved  from  the  center  position  to  any  of  the  other  LED's  and  then 
returned  to  the  center  position.  The  order  of  the  target  posi¬ 
tions  as  well  as  the  time  interval  between  displacement  was 
randomized.  Data  were  only  recorded  for  the  initial  displacement 
from  the  center  position. 

Horizontal  eye  movements  were  recorded  from  each  eye  using  an 
infrared  reflection  signal  from  the  anterior  surface  of  the 
cornea-scleral  interface  with  instrumentation  described  by 
Engelken  et  al.  (18)  Signals  for  both  eyes  tracking  were  digi¬ 
tized  using  the  analog/digital  converter  of  a  DECLAB  PDP  11/34 
computer  and  stored  in  disk  memory.  These  signals  were  sampled 
at  a  rate  of  1000  samples  per  second  for  one-half  second. 


All  parameters  in  the  model  are  estimated  with  a  conjugate 
gradient  search  program  which  minimizes  the  integral  of  the 
absolute  value  of  the  error  squared  between  the  frequency 
response  of  the  model  and  the  data.  (19)  Additionally,  the 
latent  period,  the  time  interval  between  the  target  movement  and 
the  start  of  the  eye  movement,  is  also  estimated.  A  low-pass 
filtered  pulse  step  signal  is  assumed  for  the  agonist  and  anta¬ 
gonist  muscle  forces.  (20)  Both  muscles  has  the  same  activation 
and  deactivation  time  constants.  Physiological  data  is  used  to 
estimate  all  of  the  initial  parameter  estimates.  (12)  (21) 

The  estimated  eigenvalues  of  the  fourth-order  oculomotor  plant 
correlated  well  with  their  physiologically  derived  values.  Table 
I  shows  the  estimated  parameter  values  describing  both  the 
fourth-order  and  the  second-order  oculomotor  plant  model. 

.  TABLE  I 

Parameter  values  for  the  second-order  and  the  fourth-order 
oculomotor  plant  models 

PQ  =  7.56  x  108 
P1  =  5.152  x  107 
P2  =  5.948  x  105 
P3  =  1.726  x  103 
dA  =  4.890  x  107 
d2  =  4.213  x  105 

Equations  6  and  7  are  used  to  compute  the  parameter  values  for 
the  second-order  model.  Close  agreement  between  the  predicted 


saccadic  eye  movement  response  for  both  the  fourth-order  and  the 
second-order  models  and  the  data  is  seen  for  input  target  dis¬ 
placements  from  -20  to  20  degrees.  Figure  2  illustrates  the 
observed  saccadic  eye  movement  response  for  a  10  degree  movement 
and  the  predicted  eye  movement  responses. 

• 

4.  DISCUSSION 

The  objective  of  this  study  was  the  development  of  a  physiolog¬ 
ically  based  second-order  model  of  the  saccadic  eye  movement 
system  that  is  potentially  of  value  in  evaluating  normal  and 
pathological  oculomotor  function.  We  feel  we  have  accomplished 
this  objective  as  illustrated  in  the  previous  sections.  The 
model  is  straightforward  without  obscuring  details  and  utilizes  a 
minimum  of  assumptions  concerning  the  muscle  forces.  The  second- 
order  model  performed  as  well  as  the  more  complex  fourth-order 
model  and  closely  matches  the  eye  movement  data. 

The  derivation  of  the  model  is  sufficiently  general  that  no 
assumptions  are  made  for  the  saccadic  muscle  forces ,  and  the 
model  is  simple  enough  to  be  easily  incorporated  in  oculomotor 
studies  with  little  loss  of  precision.  For  example,  investi¬ 
gators  may  utilize  this  model  when  examining  alternatives  to  the 
pulse-step  muscle  force,  a  muscle  force  that  is  a  crude  approxi¬ 
mation  to  the  actual  waveform.  (14)  The  mechanical  components  of 
the  model  are  connected  such  that  the  rate  of  change  of  the 
muscle  forces  together  with  the  actual  muscle  forces  drive  the 


LEFT  EYE  RESPONSE  IN  DEGREES 


extremely  high  initial  burst  of  motoneuronal  firing  causes  the 
rate  of  change  muscle  forces  to  move  the  eyeball  much  more 
quickly.  Additionally,  the  rate  of  change  of  the  muscle  forces 
act  as  a  dynamic  brake  which  is  slowly  reduced  to  zero  as  the 
eyeball  reaches  its  destination. 

Both  the  second-order  model  and  the  fourth-order  model  are 
derived  such  that  the  initial  conditions  for  the  series  elastic 
and  length  tension  elements  are  included  in  the  estimates  of  the 
initial  agonist  and  antagonist  muscle  forces  and  eye  position. 
As  a  result,  all  parameters  and  initial  conditions  in  these 
models  are  either  directly  measured  or  estimated  from  physiologi¬ 
cal  data. 

While  the  second-order  model  is  a  approximation  of  the  fourth- 
order  model  no  loss  of  precision  is  evident  in  the  simulation 
studies  as  illustrated  in  Figure  1,  and  both  predictions  match 
the  data  extremely  well.  The  approximation  involved  minimizing 
the  frequency  response  difference  between  the  second-order  and 
the  fourth-order  models.  Naturally,  the  simulation  results  for 
the  second-order  and  fourth-order  model  will  diverge  for  the 
higher  order  derivative  terms  such  as  eyeball  acceleration  and 
jerk.  The  divergence  is  especially  prominent  during  the  muscle 
force  transition  period,  that  is,  during  the  initial  rise  and 
fall  of  the  pulse  signal.  This  is  due  to  the  "low  pass  filter¬ 
ing"  effect  of  the  extra  two  poles  in  the  fourth-order  model. 
One  wonders  whether  such  differences  are  a  characterization  of 


muscle  forces  that  are  not  truly  representative  of  their  wave¬ 
forms,  or  an  inadequate  description  of  the  oculomotor  plant 
function . 
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APPENDIX 

OCULOMOTOR  PLANT  DERIVATION 


A  linear  oculomotor  plant  model  consistent  with  physiologic  evi 


dence  is  to  be  developed  in  this  section.  It  is  assumed  that  the 
lateral  and  medial  rectus  muscles  are  identical,  such  that  xpi 
equals  xp4  and  Xp3  equals  xp2.  Furthermore  when  binocular 
tracking  of  targets  is  involved,  one  eye,  and  perhaps  both  eyes 
may  not  be  at  the  primary  position  at  the  onset  of  movement. 
That  case  will  also  be  considered  here.  The  following  relation¬ 
ships  evident  from  Figure  1  exist  among  variables  x^,  x^,  9  and 
x. 


*!  -  *  -  xpi ' 

CA1) 

x2  =  x  +  xp4. 

(A2) 

and 

0  =  57.296  x/r. 

(A3) 

The  initial  condition  is  assigned  to  variables  x^,  x4,  and  0 
through  x;  that  is,  variable  x  at  time  0,  x(0),  and  equations  A1 
through  A3.  Initial  conditions  are  also  assigned  to  variables  x2 
and  x3  as  x2(0)  and  x3(0).  At  primary  position  x(0)  equals  zero, 
x2  equals  ~xp2  and  x3  equals  xp3. 


Equations  1  through  3  are  written  as 


FAG  =  KSTX2 

+  EAG*2 

+  KSE(xpl 

-  x)  , 

(A4) 

fant  "  kse(x 

+  xpl) 

“  KSTX3 

BANTX3  r 

(A5) 

and 

KSE(x2  +  X3 

-  2x)  = 

Jx  +  Bx  + 

Kx , 

(A6) 

where 


KST  =  KSE  +  KLT' 

J  =  57.296  Jp/r2 

B  «  57.296  Bp/r2, 

K  <=  57.296  K_/r2. 


and 


Equations  A4  through  A6  are  written  as  one  equation  by  taking  the 
Laplace  transform  of  each  equation  and  eliminating  variables  x2 
and  x3  yielding 


KSE(SBANT  +  KST)L{FAG}  “  KSE(SBAG  +  KST)L{FANT}  +  IC  ~ 
(C4S4  +  C3S3  +  C2S2  +  CjS  +  CQ)X(S), 


(A7) 


where 


C4  "  JBANTBAG' 

C3  c  jkst(bag  +  bant^  +  bbantbag' 

C2  c  jkst  +  bkst(bag  +  BANT)  +  BAGBANT(K+2KSE)  r 
ci  =  bkst  +  (bag  +  bant}  *kkst  +  2KSEKST  "  kse  *  * 
co  =  kkst  +  2ksekstklt' 


IC  "  kse  <bag  "  bant)xpi  ~  ksebag(sbant  +  KST)x2(0)  + 
KSEBANT(SBAG  +  KST)XP3  +  (P4s3  +  P3S  +  P2S  +  BKST2)x(0) * 


(A8) 


Simplification  of  the  expression  for  IC  is  obtained  from  the 
steady  state  conditions  of  equations  A4  through  A6r  that  is 

PAG(0)  =  KSTX2(0)  +  KSE<XP1  "  x(0))'  '  (A9) 

FANT*°^  =  KSE*X*B^  +  XP1*  *”  KSTX3^0^'  (A10) 


x2(0)  +  x3(0)  =  (K/Kse  +  2)  x(0> . 


(All) 


These  three  equations  are  used  to  completely  remove  the  initial 
conditions  x2(0)  and  x3<0)  from  equation  A8,  resulting  in 


ic  =  k_(br„fr„(o>  -  bantfant(o>)  +  (C4SJ  +  C,SJ  +  C^S  +  C,  - 


SE  AG  AG 


C0(BAG  +  bant)/kst)x(0> * 


(A12) 


The  equation  describing  the  oculomotor  plant  is  now  written  as 


kse*sbant  +  kst)l{fag}  "  kse(sbag  +  kst)l{fant*  + 
kse(bagfag(0>  ~  bantfant(0>)  “  co(bag  +  bant)x<0,/kst  c 

(C4S4  +  C3S3  +  C2S2  +  C-jS  4  CqJXCS)  -  (C4S3  4  C3S2  4  C2S  4 
Cj)  x(0) .  (A13) 


Further  simplification  of  equation  A13  is  obtained  by  assuming 
that  the  eye  is  undergoing  a  horizontal  saccadic  movement. 
Regardless  of  the  type  of  force  used  to  drive  the  eyeball  to  its 
final  position,  for  example,  a  pulse-step  force,  there  is  a  con¬ 
stant  force  acting  on  the  eyeball  until  it  moves.  To  be  more 
explicit,  the  agonist  and  antagonist  forces  are  decomposed  into 


the  following  functions 


„  f  fag(0) 
l  fag'  +  1 

r  fant<o) 

t  FANT*  + 


+  fag(0> 


fant<0) 


t  <  t. 


t  *  t . 
t  <  t^ 
tit; 


(A14) 


(A15) 


where  td  is  the  time  at  which  the  eye  starts  to  move.  FAG*  and 

r  •  are  zero  for  t  less  than  t..  Substituting  expressions  A14 
ANT  c 

and  A15  into  equation  A13  yields 

K_„K  (F ( 0 )  —  F.Mrtl(0))/S  4  Kcp(SB.Mrp  4  L{F._  }  - 


*'SE  ST '  AG  ANT 

kse(sbag  +  kst>l{fant'} 


’  '  SE  ANT  '  ST - AG 

(C4S4  4  C3S3  4  C2S2  4  CjS  4  CQ) 


x(s)  -  (c4s3  4  c3sz  4  C2S  4  Cj)  X (0)  . 


(A16) 


This  result  follows  since 


SKSEBANTFAG(0)/S  “  SKSEBAGFANT(0,/S  +  KSE(BAGFAG{0) 

bantW0,)  -  co(bag  +  bant)  x(0)/kst* 


Note  that  FAG'  and  FANT'  are  identical  to  the  original  agonist 


and  antagonist  forces,  respectively,  except  that  the  initial  for¬ 
ces,  FAq(0)  an<*  PANT^  9  are  ce®°ved.  Equation  A16  is  rewritten 
by  utilizing  expressions  A14  and  A15 

KSEKSTL{PAG  "  PANT}  +  SKSEL{BANTPAg'  “  BAGFANt’}  “ 
c4s4  +  c3s3  +  c2s2  +  C^S  +  CQ)  X (S)  - 

CC4S3  +  C3S2  +  C2S  +  Cx)  x(0)  .  (A17) 

By  transforming  equation  A17  back  into  the  time  domain  and  noting 
that 


and 


dP 

ZLAS 

dt 


dt 


9 


we  have 


dF 


dt 


dP  ' 
—ANT 
dt 


9 


kse(kst(fag  pant>  +  bantpag  “  bagfant)  85 

C4x  +  C^x  +  C2x  +  CjX  +  CqX. 

Since  data  are  available  for  the  angle  of  the  eyeball  with  res¬ 
pect  to  the  primary  position,  the  angle- 6  is  used  instead  of  x 
which  yields 


S(KST(FAG  “  PANT)  +  BANTFAG  “  BAGFANT>  = 
e  +  p3g  +  p2e  +  Pj©  +  po0. 

where 


S  =  57.296  KgE/rC4, 
P3  “  C3/C4' 

P2  "  W 

P1  "  Ci/C4' 
and  P0  -  C0/C4. 


(A18) 


ESTIMATION  OF  SACCADIC  EYE  MOVEMENT  MUSCLE  FORCES 
USING  SYSTEM  IDENTIFICATION  TECHNIQUES 

John  D.  Enderle 

Department  of  Electrical  and  Electronics  Engineering 
North  Dakota  State  University 
Fargo,  North  Dakota  58105 

System  identification  techniques  were  used  to  estimate  muscle  forces 
during  saccadic  eye  movements.  Horizontal  eye  movements  were  digitally 
recorded  from  infrared  reflection  signals  off  the  anterior  surface  of  the 
cornea.  The  measured  transfer  function  of  the  oculomotor  system  was  calcu¬ 
lated  from  the  fast-eye  frequency  response  to  a  step  target  displacement. 

A  fourth-order  linear  model  described  the  mechanical  components  of  the  eye¬ 
ball  and  muscle  tissue  and  a  low-pass  filtered  pulse-step  signal  described 
the  developed  muscle  force.  Physiological  data  was  used  to  estimate  the 
initial  parameter  estimates,  including  initial  conditions  for  the  series- 
elastic  and  length-tension  elastic  elements.  Final  parameter  estimates  were 
calculated  for  the  model  by  using  a  conjugate  gradient  search  program  which 
minimizes  the  integral  of  the  absolute  value  of  the  error  squared  between 
the  model  and  the  data.  The  estimated  agonist  pulse  has  a  duration  approxi¬ 
mately  equal  to  the  length  of  the  saccade  and  a  target  displacement  dependent 
magnitude  of  less  than  100  grams  for  all  eye  movements.  Time  constants  for 
activation  and  deactivation  were  estimated  to  be  3  and  4  msec  respectively. 
Steady  state  muscle  force  estimates  appeared  to  be  in  good  agreement  with 
the  length- tension  curves. 
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ESTIMATION  OP  SACCADIC  EYE  MOVEMENT  MUSCLE  FORCES 
USING  SYSTEM  IDENTIFICATION  TECHNIQUES 


John  D.  Enderle 

Department  of  Electrical  and  Electronics  Engineering 
North  Dakota  State  University 
Fargo,  North  Dakota  58105 


INTRODUCTION 

The  present  investigation  utilizes  system  identification  techniques  to  estimate 
muscle  forces  during  horizontal  saccadic  eye  movements.  Recently,  Bahill, 

Latimer,  and  Troost  (1980)  presented  a  fourth  order  linear  oculomotor  model  that 
accurately  described  saccadic  eye  movement.  Parameters  were  estimated  using  a 
time  domain  parameter  estimation  routine  (Latimer  and  Bahill,  1979).  The  large 
and  brief  estimated  agonist  muscle  force,  however,  was  not  consistent  with  phys¬ 
iological  evidence  (Robinson,  1981).  This  project  uses  frequency  response  methods 
to  estimate  muscle  forces  on  a  fourth  order  oculomotor  plant  during  saccadic  eye 
movements. 


METHODS 

The  subject  was  seated  before  a  target  display  of  nine  small  red  light-emitting 
diodes  (LED)  each  separated  by  five  degrees.  The  subject  was  instructed  to  follow 
the  jumping  target  which  moved  from  center  position  to  any  one  of  the  lED's  and 
then  returned  to  center  position.  The  order  of  the  target  positions  as  well  as 
the  time  interval  between  displacements  was  random.  Data  was  recorded  only  for 
the  initial  displacement  from  center  position.  Horizontal  eye  movements  were 
recorded  from  the  infrared  reflection  signal  off  the  anterior  surface  of  the 
cornea  from  instrumentation  similar  to  that  described  by  Young  and  Sheena  (1975) . 
Signals  for  both  eyes  tracking  were  recorded  using  the  analog/digital  converter 
of  a  DECLAB  PDP  11/34  computer  and  stored  in  disk  memory.  These  signals  were 
sampled  at  a  rate  of  1000  samples  per  second  for  one-half  second. 


Model  of  the  Oculomotor  Plant 

The  diagram  in  Fig.  1  illustrates  the  mechanical  components  of  the  oculomotor 
plant  for  horizontal  eye  movements.  The  lateral  and  medial  rectus  muscles  were 
modeled  as  a  voigt  element  in  series  with  spring  KgE.  The  passive  elasticity  from 
each  muscle  was  included  in  spring  Kp.  This  model  is  consistent  with  the  linear 
homeomorphic  model  presented  by  Bahill,  Latimer,  and  Troost  (1980)  except  for  the 
treatment  of  the  initial  conditions.  The  muscles  are  shown  to  be  extended  from 
equilibrium  position  at  primary  position  consistent  with  physiological  evidence 
(Robinson  and  co-workers,  1969).  The  following  fourth  order  model  was  derived 
for  the  oculomotor  system  shown  in  Fig.  1. 


Pig.  1.  Mechanical  components  for  oculomotor  plant. 


kse  <banfag  ~  bagfan  +  (kse  +  klt)  (f^  “  fan))  ■  p*x  +  p3x  +  p2x  +  pxx  +  p0x,  (l) 

where  Pi»  «=  D  JpBanBAG» 

P9  s  D  Jp  (Ban  +  bAG^  +  B  BpBanBag t 

P2  -  D  JpKST2  +  D  BpKgx  (Bam  +  Bag)  +  BANBAg  (O  KP  +  2KSe)  , 

Pi  -  D  BpKgT2  +  (D  KPKst  +  2KSEKST  -  KSE2)  (B^  +  Bag)  ' 

Pp  =  D  KpKST2  +  2KSEKSTKLT, 

*ST  ”  kSE  + 

r  =  radius  of  the  eyeball, 
and  D  »  57.296/r2. 

Equation  1  was  derived  such  that  the  initial  conditions  for  the  series-elastic 
and  length-tension  elastic  elements  were  included  in  the  initial  agonist  and  anta¬ 
gonist  forces  and  eye  position.  This  result  means  that  all  parameters  and  initial 
conditions  can  be  estimated  from  physiological  data. 


System  identification  Technique 

The  system  identification  technique  is  a  frequency  response  method.  Ideally, 
the  transfer  function  is  equal  to  the  ratio  of  the  Fourier  transform  of  the  out¬ 
put  to  the  system  input.  For  the  oculomotor  system,  the  transfer  function  was 
calculated  from  the  fast  eye  response  to  a  step  target  displacement.  Because  of 
unequal  time  delays  and  eye  displacements  from  saccade  to  saccade  for  the  same 
target  displacement,  it  was  not  possible  to  use  averaging  techniques  to  reduce 
the  effects  of  measurement  noise.  Fortunately,  the  measurement  noise  was  small 
relative  to  the  input  and  output  signals,  so  the  second  method  described  by  Otnes 
and  Enochson  (1972)  was  used.  The  fast  eye  response  measurements  were  first 
filtered  using  a  Butterworth  low-pass  filter  with  a  half-power  point  at  125  hertz. 
Transforming  the  filtered  measurements  directly  by  the  Fast  Fourier  algorithm 
resulted  in  distortion  due  to  truncation,  since  the  signal  did  not  go  to  zero  at 
steady  state.  This  was  circumvented  by  subtracting  the  steady  state  value  from 
each  sample,  passing  this  signal  through  a  Kaiser  window,  packing  with  zero's, 
and  Fast  Fourier  transforming  the  modified  signal.  The  Fourier  transform  of  the 
fast  eye  response  was. equal  to  the  Fast  Fourier  transform  of  the  modified  signal 
plus  the  Fourier  transform  of  the  unit  step  with  amplitude  equal  to  the  steady 
state  value.  The  input  signal  was  the  Fourier  transform  of  the  unit  step  function 
with  amplitude  equal  to  the  steady  state  value.  Parameter  estimates  for  the  oculo¬ 
motor  model  were  calculated  using  the  conjugate  gradient  search  program  similar 
to  Seidel's  (1975),  which  minimizes  the  integral  of  the  absolute  value  of  the 
error  squared  between  the  model  and  the  data.  A  low-pass  filtered  pulse-step 
signal  described  the  developed  muscle  force.  Physiological  data  was  used  to 
estimate  the  initial  parameter  estimates  (Collins,  O'Meara,  and  Scott,  1975} 

Latimer  and  Bahill,  1979;  Robinson  and  co-workers,  1969). 


RESULTS 

A  close  agreement  between  the  predicted  saccadic  eye  response  and  the  data  was  seen 
for  input  target  displacements  from  -20  to  20  degrees.  All  parameters  in  equation 
were  estimated  with  the  parameter  estimation  routine .  The  eigenvalues  of  the 
oculomotor  plant  correlated  well  with  their  physiological  data  derived  values. 

The  estimated  agonist  pulse  had  an  effective  duration  equal  to  the  length  of  the 
saccade.  The  actual  target  displacement  dependent  pulse  width  was  less  than  the 
saccade  duration,  but  because  of  the  low  pass  filtering,  effectively  equalled  the 
saccade  duration.  The  agonist  pulse  magnitude  was  also  dependent  on  target  dis¬ 
placement,  with  estimated  values  of  less  than  100  grams  for  all  eye  movements. 

The  antagonist  muscle  force  was  assumed  to  have  the  same  pulse  width  as  the 
agonist  muscle  force  and  a  magnitude  of  zero.  Steady  state  or  step  muscle  force 
estimates  appeared  to  be  in  good  agreement  with  the  length- tens ion  curves  reported 
by  Robinson  (1981).  Time  constants  for  activation  and  deactivation  were  estimated 
to  be  approximately  10  and  5  msec,  respectively.  Shown  in  Fig.  2  and  3  are  the 
time  and  frequency  response  for  a  10"  movement. 


DISCUSSION 

System -identification  techniques  were  used  to  estimate  the  eigenvalues  of  the 
oculomotor  plant  and  muscle  forces  during  saccadic  eye  movements.  While  the  eigen¬ 
values,  pulse  width,  pulse  height,  and  steady  state  muscle  forces  appear  to  be 
consistent  with  the  physiological  data  and  estimates  of  other  investigators  (Bahill 
Latimer,  and  Troost,  1980;  Collins,  O'Meara,  and  Scott,  1975;  Robinson  and 
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Pig.  2.  Time  response  for  a  10°  eye  Movement. 


co-workers,  1969),  the  estimates  of  the  time  constants  appear  inconsistent  with 
the  observations  reported  by  Robinson  (1981).  One  explanation  is  that  the 
developed  muscle  force  is  not  a  low-pass  filtered  pulse  step  but  a  more  complex 
signal.  While  the  time  response  indicates  an  excellent  fit  between  model  and 
data,  there  appears  a  more  significant  lack  of  fit  in  the  magnitude  frequency 
response.  This  lack  of  fit  may  be  indicative  of  a  more  complex  muscle  force 
signal. 
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Fig.  3.  Frequency  response  for  a  10°  eye  movement. 
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